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ABSTRACT 

In this work we build a detailed dynamic model for a SO galaxy possibly hosting 
a central massive dark object (MDO). We show that the photometric profiles and the 
kinematics along the major and minor axes, including the I13 and hi profiles, imply the 
presence of a central MDO of mass M M do ~ 1-5 - 9.7 • 1O 8 M , i.e. 0.3 - 2.8% of the 
mass derived for the stellar spheroidal component. Models without MDO are unable 
to reproduce the kinematic properties of the inner stars and of the rapidly rotating 
nuclear gas. 

The stellar population comprise of an exponential disc (27% of the light) and 
a diffuse spheroidal component (73% of the light) that cannot be represented by a 
simple dc Vaucouleurs profile at any radius. The M/L ratios we found for the stellar 
components (respectively 3.3 and 6.6) are typical of those of disc and elliptical galaxies. 

Key words: galaxies: individual: NGC 4350 - galaxies: kinematics and dynamics - 
galaxies: structure - galaxies: nuclei - dark matter 



1 INTRODUCTION 

There is increasing evidence that most galaxies with a large 
spheroidal component host a Massive Dark Object (MDO) 
at their centre, wi th masses ranging from ~ 10 8 Mq to 2 ■ 
10 10 M o ( |Ho, 199§ ). 

T his ev idence was first reviewed by Kormendy & Rich- 
stone (1995), suggesting that in kinematically hot galaxies 
(ellipticals and SOs) the mass of the central object A/mdo 
correlates with the mass of the hot stellar component M jp u 
For the 



On the other hand, using several hundred rotation 
curves, Salucci et al. (200C) recently put stringent upper lim- 
its on the mass of the MDOs resident in late-type galaxies, 
showing that they must be on average at least 10-100 times 
smaller than the ones claimed for ellipticals. If one assumes 
that MDOs are the remnants of past QSO activity, this im- 
plies that the contribution to the shining phase of quasars 
by BHs hosted in late-type spirals is totally negligible. 

Finally, recent r esults ( Ferrarese fc Merritt, 200o| ; |Gcl:> 



ratio x = log( K-hmr> /A'f CFh ) they found a. Oa.ns- 



hardt et al, 2000a) showed that the central MDO mass 



sian distribution with average value x = —2.5 and variance 
a — 0.2. More recently, assuming iso tropy in the veloc- 
ity dispersion tensor, Magorrian et al. (1998) exploited the 



high resolution of HST photometry and ground based spec- 
troscopy to estimate the MDO mass of 36 E/SOs, finding a 
similar correlation x = —2.28 ± 0.5. 

Exploiting the data from a sample of 30 galaxies, for 
which at least two inde pende nt MDO mass estimates are 
available, Salucci et al. (1999) found an average valu e x = 
—2.60 ±0.3, quite close to that obtained by Ho et al. (199!:) 
and Ford et al. (1999). This result has been used to de- 



termine the mass function of MDOs hosted in spheroidal 
galaxies and turns out to be fully consistent with the mass 
function of accreted matter required to fuel QSO activity. 
This supports the idea that MDOs are the remnant BHs of 
a past era of nuclear galactic activity. 
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strongly correlates with the galax y velocity dispersion. In 
particular Gebhardt et al. (2000a) pointed out that ellip- 
ticals and bulges fall again in a fundamental plane when 
they are located in a four- dimension plane with coordinates 
(Mbh — logL — log a, log R e ). This property implies some 
fundamental connection between the black hole and the 
bulge. 

In this context it is interesting to investigate the role 
of SO and early-type disc galaxies. In view of their mas- 
sive bulge, they must be considered as a possible location 
for QSO remnants. However, the search for central massive 
objects in SO is hampered, as: (1) the bulge/disc photomet- 
ric decomposition is often uncertain; (2) these objects have 
weak Ha lines and, consequently, the rotation curve inside 
the central 0.5 kpc cannot be accurately determined; (3) the 
rotation curve is too disturbed by the strong asymmetric 
drift effects to be representative of the circular velocity. 

Detection of MDOs in this kind of galaxy requires kine- 
matic data of the innermost regions, and a detailed dynamic 
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Table 1. Properties of NGC 4350 ( gully, 1988[ ). (1) Hubble type 
T; (2) absolute blue magnitude \ (3) total blue luminosity Lg; 
(4) position angle P.A.; (5) observed apparent axial ratio at D25; 
(6) distance D; (7) conversion scale from arcsec to parsec; (8) 
observed diameter at 25 mag/arcsec 2 isophote in blue. 
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dispersion profile suggest the presence of a central concen- 
tration of mass. 

Loyer et al. (1998) modeled this galaxy in some detail 
using a Multi- Gaussian Expansion method. In this paper, 
the photometric and kinematic data are reproduced with a 
different method and the subsequent analysis is extended to 
include the gas rotational velocity and the profiles of the 
/13, /14 paramete rs. As pointed out by several authors ( Van 



Der Marel, 1999 



Cretton & Van Den Bosch, 1999 



Gebhardt 



[T] [mag] [Mq] [°] 



[Mpc] [pc///] H_ et < 20001)[ ). the use of the higher moments of the line-of- 
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sight velocity distribution is very relevant in determining the 
orbital anisotropics, thus reducing the uncertainties on the 
MDO mass. 



model. This approach has proved effective in detecting mas- 
sive central dark objects in a few S O/Sa g alaxies, including 



NGC 4594 ( [Kormendy et al, 1996f), M31 ( [Tsvetanov et al. . 
199cOJGC 3115 (lEmsellem et al, 1999| ), and NGC 4342 
(Cretton & Van Den Bosch, 199£). Given the highly com- 



pact stellar structure, central MDOs can only be detected 
through accurate reconstruction of the galactic components, 
including the bulge/disc decomposition. 

The goal of this work is to investigate the possible pres- 
ence of an MDO in the nucleus of the SO galaxy NGC 4350, 
by ap plying the dynamical model of Pignatelli & Galletta 
(1999) to recent data. This work has been set out as follows: 
in Section ^ we present the available data; in Section ^ we 
give a brief summary of the hypothesis and characteristics 
of the model; in Section ^ we derive the best-fit parameters 
of the model; Section. [B| is devoted to the discussion of the 
results and conclusions. 

In this paper we assume a Hubble constant Ho = 75 km 
s" 1 Mpc" 1 . 



2 THE DATA 

NGC 4350 is classified as a classical SO galaxy in the Nearby 
galaxy catalog (Tully 1988; see Table 0). It is observed almost 
edge-on (from its axial ratio one can infer an inclination 
angle i > 62°). 

For the photometry, we adopt the R-band surface 



brightness profiles published by Seifert and Scorza (1996). 
According to these authors, NGC 4350 can be decomposed 
into disc and bulge components. The disc is probably expo- 
nential, significantly decreasing beyond 40 arcsec; the bulge 
on the other hand dominates even in the outer parts of the 
galaxy, as indicated by the substantial drop in the elliptic- 
ity profile. No ring-like component, or inner stellar disc has 
been detected. 

For the kinematics, we adopt Fisher's detailed stellar 
kinematical o bservations along the major and minor axis 
(Fisher, 1997). These include the rotational velocity, the 
velocity dispersion profile, and the pro files of the Gauss- 
Hermite higher-order moments /13, /14 (Van Der Marel & 



Franx, 1993). Moreover, they also include emission-line gas 



kinematics of a rapidly rotating nuclear disc. The total num- 
ber of indipendent measurements to be simultaneously fitted 
by a model is 280. 

The high rotational velocity of the nuclear gaseous disc, 
the rapid increase of the stellar rotation curve in the first 2" 
(~ 160 pc) and the corresponding high peak in the velocity 



3 THE GALAXY MODEL 

In order to investigate the gas and stellar kinematics we 
use the se lf-con sistent dynamical technique of Pignatelli and 
Galletta ( 1999| ). We give here a brief summary of the tech- 
nique and the general assumptions made. 

A galaxy can be described by superposition of different 
components. For each component, we separately assume: 

- the density distribution is oblate, without triaxial struc- 
tures; 

- the isodensity surfaces are similar concentric spheroids; 

- the surface density profile follows a simple r 1 ^ 4 or an 
exponential law; 

- the local velocity distribution is Gaussian; 

- the velocity dispersion is isotropic (ay = oe — <r z ); 

- the mass-luminosity ratio is constant with radius; 

Our model does not consider the possible presence of triaxial 
structures (bar; triaxial bulge; tilted component; warp) or of 
anisotropy in the velocity distribution (with predominance 
of radial or tangential orbits). 

The model has 4n + 2 free parameters, where n is 
the number of adopted components: namely the lumi- 
nosity, scale length, mass-luminosity ratio and flattening 
{Ltot, r e , M/L, b/a} of each component plus the inclina- 
tion angle of the galaxy and the mass of the central massive 
dark object Mm do- In principle, photometry can be used 
to constrain all these parameters except the M/L ratios and 
the value of Mm do, which must be derived through the 
kinematics. 



4 COMPARISON BETWEEN MODEL AND 
DATA 

4.1 Photometry 



Seifert and Scorza (1996) have performed detailed CCD pho- 
tometry in the R-band for a sample of 16 early-type disc 
galaxies, including NGC 4350. For each galaxy, they de- 
termined the radial profiles for surface brightness and el- 
lipticity, together with a bulge/disc decomposition without 
special assumptions on the radial density profiles of the two 
components. According to the authors, the disc of NGC 4350 
extends up to the innermost regions of the galaxy and follows 
approximately an exponential law; the spherical component 
dominates both in the nucleus and in the outer parts, as 
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Figure 1. Photometric decomposition of NGC 4350. The parameters of the different components are shown in Table Data are taken 
from Seifert and Scorza (L996). The dotted lines represent the separate contributions of the three different components to the global 
luminosity profile. 



scale radius 



axial ratio 



Luminosity 



r& r d v h 

2±0'.'l 15±0'.'2 40 ±( 



(b/a) 6 (b/a) d (b/a) h 

l.±0.03 0.1 ±0.01 0.95 ±0.06 



Lb L d L h 

24.5 ±4% 26.5 ±1.5% 49 ± 5% 



85±5° 



Table 2. Parameters of the decomposition shown in Fig. [l|. The errors define the la level in the parameter space on the basis of a 
reduced x-square analysis. 



indicated by the strong external decrease in the ellipticity 
profile (see Fig. [j]). 

In order to describe the photometry of this galaxy we 
have adopted a three-component model: 

(1) an almost edge-on exponential disc component; 

(2) a diffuse spheroidal component, with a deVaucouleur 
surface brightness profile, extending well beyond the disc 
itself, identified hereafter as "stellar halo"; 

(3) a second spheroidal component, very compact and 
with an exponential profile in the inner 2" (see Fig. hi), which 
we identify as "bulge". 



photometric data are presented in Table |^, together with 
their statistical uncertainties. 

The rapid increase in the ellipticity from very small val- 
ues to a maximum of e ~ 0.7 at r ~ 30" (see Figure [l]) 
and its slow decline in the outer parts can be explained if 
at least three components are used. The two "spheroidal" 
components may represent two physically different stellar 
systems, or, since the best-fit Battenings are very similar, 
a single, almost spherical component with a density profile 
not represe nted by the de Vaucoleur's law (see discussion in 
the Section 5.2). 



The parameters derived by performing a best fit to the 
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Figure 2. Data from Fisher (1997) for the major axis, compared with our best-fit models with (solid line) and without (long-dashed 
line) a central MDO. The squares and triangles represent data respectively from the approaching and the receding side. From the top: (a) 
stellar (filled symbols) and gas (open symbols) rotational velocities. The latter is also compared with the circular velocity inferred from 
the two models; (b) velocity dispersion profile; (c) /13 radial profile; (d) Z14 radial profile. Data beyond 2.5 kpc are not shown because 
they were not included in the fit. 



4.2 Kinematics 



Fisher ( 1997 ) has published ground-based kinematics of 18 
early-type galaxies, including NGC 4350. The data includes 
the stellar rotation curve, velocity dispersion profile and the 
line-of-sight velocity distribution described in terms of the 



Gauss-Hcrmitc expan sion series terms /13 and /14 (Van Der 
Marel & Franx, 1993). Emission line spectroscopy reveals a 



disc of rapidly rotating gas within the inner 3", which ap- 
pears kinematically decoupled from the stellar component. 
In fact, while the gas is likely to follow circular orbits, the 
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Figure 3. Data from Fisher (1997) for the minor axis, compared with our best-fit model with (solid line) and without (long-dashed line) 
a central massive dark object. The notation used for the symbols is the same as in Fig. |^. 

Most parameters of the model have been fixed by fitting 
the photometry; the kinematic data are fitted by varying 
only the M/L ratios of the different stellar components and 
the value of the MDO mass. 

A robust upper limit to the central mass of the MDO 
can be obtained by assuming that the gas follows a Keplerian 
rotation curve. After convolution with seeing and slit effects, 
we obtain Mmdo < 1-5 • 1O 9 M0. However, since part of the 
gravitational force is due to the diffuse stellar component, 



stellar velocity is obtained by taking into account the asym- 
metric drift effect as 

K-ot = 'Kcirc — AVasymdrift , (1) 



where the /S.V^ syraAliit ( Binney fc Tremaine, 1987 ) is essen- 
tially a function of the velocity dispersion and is only slightly 
dependent_on the stellar distribution. The data are shown 
in Figs. 
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we allowed the value of the dark mass concentration to vary 
between zero and M max = 1.5 • 10 9 Mq. The value of the 
M/L ratios for the spheroidal components was allowed to 
vary from 1 to 10 times the corresponding value for the disc. 

We integrated t he Je ans equations as detailed in Pig- 
natelli and Galletta (1999) to reproduce the kinematic data 
(see Figs. ^-^). The best-fit stellar M/L ratios spheroidal 
components and Mm do are given in Table |§| together with 
their statistical uncertainties. 

For the MDO model, the error bars were obtained by 
computing the reduced x 2 f° r the whole grid of models, and 
rejecting the models with a reduced % 2 greater than the 95% 
confidence level. For the non-MDO model, since the reduced 
X 2 was very high for all of the models, we kept the models 
with a % 2 greater than three times the minimum value. 

From Figure |2| it is apparent that within 2.5 kpc the 
galaxy can be divided into 3 different regions: 

(1) an inner region (r < 120 pc) which features a sharp in- 
crease in both the stellar and gas rotation curve. This implies 
that this region is dynamically dominated by the presence of 
the central dark massive object, which also determines the 
central cusp in the velocity dispersion profile. 

(2) an intermediate region (120 pc < r < 700 pc) which is 
bulge-dominated. In this dynamically hot region the asym- 
metric drift effect dominates the stellar rotation curve, as 
shown by the large difference (~ 100 km/sec, i.e. a factor of 
more than 2) between the gas and the stellar velocities. The 
stellar rotational velocity reaches a plateau at about 50% of 
the circular velocity. 

(3) an outer region (r > 700 pc) which is disc dominated. 
Here the stellar rotational velocity rises again, while the 
asymmetric drift effect plays a minor role. There is a proba- 
ble flex point at r w 900 pc in the stellar velocity dispersion 
profile, as expected from the superposition of bulge and disc 
populations with different kinematical behaviour but similar 
masses. 

Remarkably, these kinematic features match exactly 
those expected from photometric decomposition. This is not 
trivial: it implies that inside 2 kpc the light traces the mass, 
thus ruling out the possibility of a sizable presence of dark 
matter. 

It is worthwhile to discuss in detail the differences be- 
tween the best-fit MDO and non-MDO models in Fig. [|||. 

First, under the assumption of M/L ratios constant with 
radius for each component, the M/L values found in the 
two cases are similar (see Table |^) within the error bars 
given by the fit. This reflects the fact that in both cases the 
kinematics beyond r — 500 pc severely constrains all the 
parameters of the luminous matter in both cases. 

The different panels of Fig. | show that the presence of 
the central massive object affects the different kinematical 
profiles up to different radii: 

• The velocity dispersion profile is affected inside a 
"sphere of influence" of radius tbh = GMbh / 'o~q. This is 
the radius at which the black hole is expected to increase 
the velocity dispersion by a factor y2- For the observed 
value of (To and the inferred value of Mbh shown in Table ^[ 
we expect reH to be about 70 pc, which is similar to what 
we see in Fig. 2b. 

• While the velocity dispersion depends on the whole 




Figure 4. Luminous mass distribution of NGC 4350, according 
to Table We plot the contribution of the inner exponential 
bulge [dotted line), the exponential disc {short-dashed line) and 
the outer deVaucouleur stellar halo (long- dashed line). The total 
mass is also shown as a solid line. Note that at r = the mass is 
greater than zero because of the contribution from the MDO. 



mass distribution, the circular velocity is only sensitive to 
the inner mass. The "radius of influence" in this case is the 
radius at which the stellar mass equals the MDO mass, and 
corresponds to about 200 pc, as can be seen in Fig. ^| At 
this radius the circular velocities of the MDO and non-MDO 
model differ by a factor y/2, as shown in Fig. 2a. Beyond 
r = 300 pc part of these differences is due to the different 
stellar M/L ratios (see Table ^). 

• The differences between the two models in the circular 
velocity are strongly amplified in the stellar rotation curve. 
For instance, at r = 300 pc, the MDO and non-MDO mod- 
els have respective circular velocities of about 200 and 180 
km/sec. Their asymmetric drift corrections AVasymdrift ac- 
count respectively for about 185 and 178 km/sec. If the ro- 
tational velocities are computed using Eq. (hh, these small 
differences are magnified as shown in Fig. |i[ 

• Finally, the MDO affects the /13 moment of the line-of- 
sight velocity distribution up to a similar radius (r w 300 
pc). The /13 signature, together with the gas rotation curve, 
are of great help in constraining the orbital anisotropies be- 
cause they rule out strongly radial orbits in the inner part 
of the galaxy. 



5 DISCUSSION 

The photometric and kinematic data along both the ma- 
jor and minor axes have been explained through the super- 
position of different components: an inner MDO of mass 
4.9 — 9.7- 10 s Mq, and three stellar components described in 
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Bulge 

Mass (1.3 ± 0.2) ■ 10 10 M G 
M/L 6.6 ±0.9 



Model with MDO 
Disc 



Halo 



(7.1 ± 0.7) ■ 10 9 M (2.7 ± 0.3) ■ 10 10 M c . 



MDO 
(7.3 ± 2.4) ■ 10 s M e 



3.3 ±0.3 



6.6 ±0.7 



Bulge 

Mass (1.5 ± 0.2) ■ 10 10 M 
M/L 7.1 ±1.0 



O 



Model without MDO 
Disc Halo 
(6.6 ±0.7) ■ 10 9 M (2.8 ±0.3) ■ 10 10 M 
3.1 ±0.3 6.8 ±0.6 



MDO 



Table 3. Masses of the different components for the best-fit models both with and without a central MDO. The error bars define the 
95% confidence region in the parameter space, under the assumption that each component has a M/L ratio constant with radius. The 
mass profile derived for the model with a central MDO is shown in Fig 



Table ^, with radial mass profiles M{r) illustrated in Fig. ^ 
Each component is discussed separately hereafter. 



5.1 The Massive Dark Object 

The rotational velocities of both gas and stars are known 
to be much more efficient in constraining the mass of the 
central MDO than the velocity dispersion. Models with- 
out a central black hole fail to reproduce the shape of the 
inner rotation curve (see dot-dashed line in Fig. |^), un- 
less unreasonable values of the bulge mass-to-light ratio are 
assumed which, however, would imply velocity dispersions 
much larger than those observed. In the inner kpc, the re- 
duced chi-square is \ ~ 1.77 for the best-fit model with 
MDO, and \ ~ W0 for the model without MDO. 

On the other hand, the errors given in Table ^| are just 
statistical errors. We must consider also possible systematic 
errors related to the assumptions used in the model. 

The main concern is related to different hypotheses on 
the phase-space distribution function of the stellar compo- 
nent. So far, the model we presented is axisymmetric, with 
an isotropic velocity distribution. Such a model is associ- 
ated with a distribution function which only depends on 
the two classical integrals of motion, i.e. energy and vertical 
angular momentum: / = f(E,L z ). However, strong radial 
anisotropy in the centre of the galaxy results in a high veloc- 
ity dispersion s, mimicking the presence of a MDO ( jGebhardt 



et al, 2000b) 



For this specific galaxy, significant anisotropies in the 
velocity distribution would certainly raise the central ve- 
locity dispersion of the non-MDO model up to the observed 
value, but: (i) the gas rotational velocity, which is unaffected 
by the isotropic hypothesis, would be too low with respect to 
the observations; (ii) the stellar rotational velocity would de- 
crease, moving away from the observations; (iii) The radial 
orbits in the central bin wound increase the hi parameter, 
in contrast with that suggested by the data. Thus, in this 
galaxy the hypothesis of radial orbits can not reduce the 
estimated MDO mass. 

In principle, we cannot exclude the possibility that, in- 
side the innermost observed kinematical point, the M/L 
ratio of the stellar component is rising sharply enough to 
mimic the presence of a massive dark object. However, the 
innermost observed kinematical point is at .Rmncr ~ O'.'G = 
66 pc, and the data requires a dynamical mass inside this 
point of 1.1-1O 9 M0, which, compared with the observed pho- 



tometry, implies an M/L ~ 28 in this region, much highe r 
than the value M/L « 11.5 found by Loyer et al. ( 1998 ). 
Such a high mass-to-light rat io seems unlikel y for any given 
choice of stellar population (Wbrthey, 1994). Adopting an 
extreme M/L « 12 for the innermost region would only re- 
duce the best-fit Mm do to 5 • 1O 8 M0, and the lower limit 
to 2.5 ■ 10 8 M Q . 

In addition, we have not forced the whole galaxy to have 
the same stellar mass-to-light ratio. Instead, we allowed dif- 
ferent M/L ratios for different components, although the 
best model fit suggests similar values for the bulge and 
stellar halo. This seems a reasonable hypothesis because it 
links every dynamical component to a given star formation 
history. The M/L ratio of our global stellar component in- 
creases with decreasing radius. 

Finally, the hypothesis of constant flattening for each 
photometric component is also a possible source of system- 
atic errors. Because of this assumption, we obtained a pho- 
tometric decomposition with small statistical errors (see Ta- 
ble ^), despite the large seeing (about 2'.'5) of the ground- 
based observations. If we relax this hypothesis, then it is pos- 
sible that the stellar component becomes flatter inside the 
seeing radius, contributing to the rotation shown in Fig. ^ 
and reducing the need for a central MDO. 

We can not reproduce this kind of behaviour in a self- 
consistent way with our model. We can, however, give an 
estimate of the systematic errors involved, evaluating a set 
of models with different ellipticities for the bulge component, 
and trying to reproduce the stellar V and a profile only in 
the inner 2". In order to fit the inner rotation curve without 
MDO, all the mass inside 2" must flattens to £ ^ 0.3. On 
the other hand with e = 0.3, the central velocity dispersion 
falls below the predictions for the no MDO case with e = 
and does not match the observations (cfr. Fig. []). Under the 
asssumptions of M/L = 12 for the bulge and e ~ 0.2 within 
2.5" , both the inner rotation velocity and dispersion velocity 
can be fitted within the error bars, if a central dark object 
with Mm do ~ 1.5 • 10 8 is present. 

Taking into account all the above listed source of un- 
certainties, the dark object mass falls in the range 1.5 <; 
M M do < 9.7 x 10 s , implying -2 < log(M MD o/M hulgc ) < 
— 1.1. This wide range is substantially due to the poor reso- 
lution of the available data. 

We can now locate NGC 4350 in the a e — Mmdo 
plane. F ollowing the same prescription as in Gebhardt et al. 
( fcOOOal ), we computed for NGC 4350 a luminosity-weighted 
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<j e = 180 km /s insi de R e . The relationship found by Geb- 
hardt et al. (2000a) would predict for this value 6.5 x 10 7 < 
Mm do < 9.7x10'. Although our best estimate is a factor of 
10 larger than the value expected on the basis of a e — Mm do 
relationship, nevertheless our lower limit is quite close to the 



We expect this difference in the M/L ratio of the two 
components to be visible in the colour distribution of the 
galaxy. There is indeed evidence for a colour gradient, going 
from (B — V) = 1.1 in the inner, bulge-dominate d region, 
to a (B — V) as 0.8 4 in the intermediate region (Prugniel 



expecte d range. In this context it is worth noticing that NGC 



Heraudeau, 1998), which qualitatively agrees with our 



4350 was selected for our analysis on the basis of the high 
velocity rotation of its nuclear gas. 

It is interest ing to compare our results with those given 
in Loyer et al. ( 1998 ). Adopting a single M/L ratio of 7 
for the whole stellar population, they found that the inner 
kinematics required a central point mass of M= 8-10 8 M Q , 
a result that agrees very well with our value of Mm do = 
7.3 ■ 10 s Mq, despite the difference in the methods adopted. 
However, the evidence for a MDO was marginal. T he dif - 
ference between our results and those of Loyer et al. ( 199£ ) 
is due to the different dataset adopted: in particular, the 
fact that we included both the gas and star kinematics, and 
especially the /13 and /14 profiles, dramatically increases the 
constraints we can put on the Mm do- This is especially 
true in the inner 5", where the shape of the line-of-sight 
velocity distribution is strongly non-gaussian (see Fig. H), 
corroborating the conclusion that the search for an MDO is 
much more effective if the analysis of the highe r orders of 
the line -of-sight veloci ty distribution is included (Cretton & 
Van Dqn Bosch, 1999|). 



5.2 The stellar component 

The photometry of this galaxy can not be described by a 
simple superposition of a deVaucouleur spheroid and an ex- 
ponential disc. The addition of a third spheroidal component 
is required. However, the fact that the mass-to-light ratios 
that reproduce the kinematical data for the two spheroidal 
components are similar suggests that there is a single, al- 
most spherical stellar population; its brightness profile fol- 
lows the deVaucouleurs law only up to the inner ~ 600 pc, 
and then starts to diverge towards a much steeper inner 
surface brightness profile. This behaviour agrees with recen t 



HST high resolution photometric data (Lauer et al., 199 



showing that many ellipticals and bulges do have inner cusps 
steeper than expected from the r 1//4 law, with a luminosity 
density profile I(r) ~ f~ 7 , wit h 7 ~ 1 for 5 pc < r < 1 
kpc. In addition, Faber et al. (1997) have shown that this 
behaviour is more evident in rapidly rotating, disky bulges 
such as the one in NGC 4350. 

A possible relation between the presence of a central 
BH and the inner surface brightness pr ofile has been ex- 



tensively investig ated by several authors (Young, 1980; Van 



Der Marel, 1999) in the framework of adiabatic black hole 



growth. The effect on the photometry, however, can only be 
seen in the "sphere of influence" of radius tbh = GMbh / 0%, 
that we showed to be about 70 pc, which is less than 1" at 
Virgo distance. Our poor angular resolution did not allow 
us to check for this effect. 

It is noteworthy that the mass-to-light ratio obtained 
for the spheroidal component ((M/L) sp h = 6.6 ± 0.7 ) is i n 
good agreement with the relation by van der Marel (1991), 



which yields (M/L) 3p h = 4.8 ± 1.5 for a galaxy of this lu- 
minosity. At the same time, the value we found for the disc 
( (M/L ) n = 3.3 ±0.3) is in agreement with Salucci & Persic 
(|l999|), that give (M/L) d ~ 2.5 ± 0.5 at this luminosity. 



results. 



5.3 The dark matter halo 

In the previous sections, we purposely limited ourselves to 
the inner regions of the galaxy (r < 2Rd), since here the 
diffuse dark matter component is clearly not required by 
the data. 

A complete and detailed investigation of the dark mat- 
ter distribution in NGC 4350 is beyond the scope of this 
paper. Nevertheless, under the hypotheses described above, 
it is possible to derive a preliminary estimate on the mass 
and distribution of dark matter needed to explain the kine- 
matical data in the outer re gions. 



The observed velocity (Fisher, 1997) is nearly circular 



and almost flat in the range 1.8 < r < 3.5 kpc. It can be 
represented by a linear law V(R) = (175 ± 5) + (3.6 ± 1.2)r , 
where r is expressed in kpc e V in km/sec. With the stellar 
distribution we derived in Table this implies that within 
the last observed kinematical point i?i as t ~ 3.5 kpc the dark 
matter amounts to M D m = 1.2 • 1O 1O M , which is 30% of 
the total mass inside this point. The DM density should 
equal that of the luminous matter at R ~ 3 kpc, so that, as 
observed in spiral galaxies (Salucci & Persic, 1999), in the 



inner regions the light traces the dynamic mass. 



6 CONCLUSIONS 

From the photometric and kinematic data available in the 
literature, we inferred the mass distribution for the SO 
galaxy NGC 4350. The data used to constrain the model 
include: R-band photometry, rotational velocity of gas and 
stars and velocity dispersion profiles along both the major 
and minor axis, and the /13 and /14 profiles, for a total num- 
ber of 280 data points. Our results can be summarized as 
follow: 



of the galaxy. A lower value for the central mass is possible 
if we allow for a very sharp increase of the anisotropy in the 
velocity distribution. However, this could be obtained only 
by neglecting the information content in the gas kinematics 
and the /13 and /14 parameters profiles. 

- This galaxy falls in the upper envelope of the Mbuige — 
Mbh or <r e — Mbh relations known so far. This is expected, 
given that it has been selected for the very fast rotation of 
its nuclear disc. 

- The estimated stellar M/L ratios (~ 6.6 for the spheroid 
and ~ 3.3 for the disc) are in agreement with the average 
values usually found for these components. In order to ex- 
plain the mass concentration in the centre of NGC 4350, 
the stellar M/L ratio has to increase in the innermost 60 pc 
from the average value of 6.6 to ~ 28. Such a high value is, 
however, ruled out by the comparison between the observed 
colours and synthetic stellar population models. 
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- Within the hypothesis of the model, the kinematical data 
between 2.5 and 3.5 kpc imply that the DM density equalizes 
that of the luminous matter at R « 3 kpc, and that the dark 
matter mass is 30% of the total within r = 3.5 kpc. 

Acknowledgements. We thank the referee, Dr. K. Geb- 
hardt, for helpful comments and suggestions that greatly 
improved the quality of the paper. 



REFERENCES 

Binncy, J. and Tremaine, 1987. Galactic Dynamics. Princeton 

University Press. 
Cretton, N. and Van Den Bosch, F. C, 1999, ApJ, 514, 704. 
Emsellem, E., Dejonghe, H., and Bacon, R., 1999, MNRAS, 303, 

495. 

Faber, S. M., Tremaine, S., Ajhar, E. A., Byun, Y. I., Dressier, 
A., Gebhardt, K., Grillmair, C, Kormendy, J., Lauer, T. R., 
and Richstonc, D., 1997, AJ, 114, 1771. 

Ferrarese, L. and Merritt, D., 2000, preprint, astro-ph/0006053. 

Fisher, D., 1997, AJ, 113, 950. 

Ford, H. C, Tsvetanov, Z., Ferrarese, L., Dressel, L., and JafTc, 
W., 1999, in: Galaxy Interactions at High and Low redshift, 
Proceedings of the IAU conference held in Kyoto, Japan, Au- 
gust 26-30, 1997. eds. D.B. Sanders and J. Barnes, Kluwcr. 

Gebhardt, K., Bender, R., Bower, G., Dressier, A., Faber, S. M., 
Filippcnko, A., Green, R., Grillmair, C, Ho, L., Kormendy, 
J., Lauer, T. R., Magorrian, J., Pinkney, J., Richstone, D., 
and Tremaine, S., 2000, preprint, astro-ph/0006289. 

Gebhardt, K., Richstone, D., Kormendy, J., Lauer, T. R., Ajhar, 
E. A., Bender, R., Dressier, A., Faber, S. M., Grillmair, C, 
Magorrian, J., and Tremaine, S., 2000, AJ, 119, 1157. 

Ho, L., 1999, in: Observational Evidence for the Black Holes in 
the Universe, p. 157. 

Kormendy, J. and Richstone, D., 1995, ARA&A, 33, 581. 

Kormendy, J., Bender, R., Ajhar, E. A., Dressier, A., Faber, S. M., 
Gebhardt, K., Grillmair, C, Lauer, T. R., Richstone, D., and 
Tremaine, S., 1996, ApJ, 473, L91. 

Lauer, T. R., Ajhar, E. A., Byun, Y. I., Dressier, A., Faber, S. M., 
Grillmair, C, Kormendy, J., Richstone, D., and Tremaine, S., 
1995, AJ, 110, 2622. 

Loyer, E., Simien, F., Michard, R., and Prugniel, P., 1998, A&A, 
334, 805. 

Magorrian, J., Tremaine, S., Richstonc, D., Bender, R., Bower, G., 
Dressier, A., Faber, S. M., Gebhardt, K., Green, R., Grillmair, 
C, Kormendy, J., and Lauer, T., 1998, AJ, 115, 2285. 

Pignatelli, E. and Galletta, G., 1999, A&A, 349, 369. 

Prugniel, P. and Hcraudeau, P., 1998, A&AS, 128, 299. 

Salucci, P. and Persic, M., 1999, MNRAS, 309, 923. 

Salucci, P., Szuszkicwicz, E., Monaco, P., and Danese, L., 1999, 
MNRAS, 307, 637. 

Salucci, P., Ratnam, C, Monaco, P., and Danese, L., 2000, MN- 
RAS, in press. 

Seifert, W. and Scorza, C, 1996, A&A, 310, 75. 

Tsvetanov, Z. I., Pei, Y. C, Ford, H. C, Kriss, G. A., and Harms, 
R. J., 1998, in: American Astronomical Society Meeting, Vol. 
193, p. 119.01. 

Tully, R. B., 1988. Nearby galaxies catalog. Cambridge and New 

York, Cambridge University Press, 1988, 221 p. 
Van Der Marel, R. P., 1991, MNRAS, 253, 710. 
Van Der Marel, R. P., 1999, AJ, 117, 744. 
Van Der Marel, R. P. and Franx, M., 1993, ApJ, 407, 525. 
Worthey, G., 1994, ApJS, 95, 107. 
Young, P., 1980, ApJ, 242, 1232. 



© 0000 RAS, MNRAS 000, 000-000 



